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TITLE OF THE INVENTION 

Method and systeGi for detecting and classiiymg the laodulation of unknown 
analog and digital telecoitmiimicati^ms sigjiEils- 

FIELD OF THE INVENTION 

The invention relates in general to tdecomsDunications and in particular to the 
classification of modulation schemes embedded in unknown telocommmications signals, 
typically: CW, FM, PSK, AM, BPSK, QPSK, 'k/A - QPSK. MPSK, NOISE, 

OTEGBRS. 

BACKGROUND OF THE IISVENTION 

AutoiDatic recognition of the modulation scheme embedded in an uxikno'wn 
received signal is an important requirement for civilian, nailitary, and goveomment 
intelligence bodies when monitoring the radio commnnication spectrum. Although the 
subject has been extensively researched for sev^al years and different approaches have 
been implemsuted or deilneated in tiieoretical papers. On? prior art has traditionally 
assujned that (a) the earner frequeney of fhe tmknovm received signal is given and has 
zero error, that (b) tiie input SignaL-to-^Noise Ratio (SNR) i$ sufficiently high to classig^ 
the modulation correctly, md that (o) the symbol transition of digitally-modulated signals 
is knovm. Furthermore, the more recent prior art approaches are limited to offline 
operation on stored signals^ Some of the rec^t metliods employ probabilistic models to 
minimize misclassification errors, which can achieve good results at SI>flas that are as low 
as 0 dB. As shown in reference [i], however, they assume a pnori fenowledge of the 
carrier phase and freqtiency, the SZSttl, and the symbol rate of the modulation, and aj^ 
often hniited to digital phase modulation schemes. 

Other approaches to autoniatic cIas$ii7cation "use statistical pattern recognition 
techniques, suoh as Artificial Neural Networks (ANN), to discenii discruninating features. 
As shown in references [2], [3], and [4], ANN d^issifiers produce reasonably good results 
under simulated conditions, but their practical behaviour is liighly dependent on the 
training sef chosen. Since they can perform learning vector quantizmion, nenral networks 
are capable of achieving an efficient class definition over a l^ge multi-dimensional 
feature space. The inherent problem however, is that ANN classifiers have difficulty 
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indeiilifying the set of meaningfiil features and to train tJie network accotdkigly. 
Furtliermore, the designer does not have much nmtpTilative ooutrol over the classification 
algoritlxni and may hay? difficulty applying a priori knowledge of the taxooiomy of tlie 
modutoicn schemes, A ns\iral netW'Ork qp^frates like a black box tbat requires a new 
training phase when new features (or signals) to be identifisd are added. ExaiTiples of 
ANN techniques are described in [4], 

Other, mors recent^ prior art approaches to mod\ilation recognition base their 
das$iftcation decisions on a momber of successive serial tests^ each yielding a binary 
output. As described in refeiCTces [2], [3], and [4], such methods give rise to a decisign 
tree in which the outcome of the first binary deoj sian forces a second binary decision 
whose outcome determines t]ie next binary dmsiou, etc. This decision tree technique 
represents an improvement over the vector-based methods described above in that it 
refines aind clarifies, in successive deoision levels, the information extracted fix>m the 
unknown input signal- Its hierarchical stnictiire allocates the computing resources .more 
efficiently. As well, the thresholds established at each decision level may be tfuickly 
modified in order to reflect operational changes. These alterations can improve 
pcrfoimance accuracy'-- Not^vithstandmg diese advantagcsj the decision tree methods 
published to-date have inherent deficiencies, namely their intolerance to eanier frequency 
enxirs, their errsLlic performance at low SlSIBs^ and their inability to cla$sify modulcition 
schemes reliably mider real-time operations. 

Other forms of classification techniques dxt described in the following US patents: 
i:i reference [5], the classifier of an IF rfgnal ta1<e5 the outputs of two separate 
demodulators (one AM, the other FM) to compute different signal statistics (or features) 
and make six binary decisions based on those statistics. It then classifies the modulating 
scheme, within a set consisting of CW, DSB, SSB, ASK, FSfC, MUX, NOISE, and 
OTHERS, by using tho whole vector of si::i binary decisions as an input to a logic ciicirit. 
The drawback to this method is that it usually performs the computation of the vecfor 
features in parallel, without any interaction between these features. It also uses a snb- 
optimimi classification circuit 

In reference [6]^ the probability distribution of the input signal amplimde is 
analyzed to estimate the mean, the variance and the aniphtude cumulative distribution. 
This information is combined with the outputs of three phase-looked 1oq5>s - one tailored 



to AM signals^ OTie to FM signals and <?ne to DSB signals- The combined. infonnatioB is 
compared witia a numb^ of tiixcsholds to form m infoirmLtion vector which is then 
compared to a pre-stored series of vectois representing the modulatioiis i^-itiiin the set CW, 
AM^ EM, DSB, SSB, PSK, The mam differetLce hetween tlnis method and the preseiit 
inventioxi is the computationally-mtensive parallel processing of the feature vector, as 
opposed to serial processing of the vector wMch is less cotBputer hungry. 

In T£^fere^oe [7], several parameters, iijcludmg the mean amplitude, the signal4o- 
aoise ratio, and the standard deviation, are computed for each of ike fteq-uenc^^ hties of ihe 
input signal's power spectrtmi, Thes^ parameters fed hi paraUel to a aeiiral network for 
the classification of the input modiilation. The main drawback of this method is that it 
performs the computation of the vector features mparaiki over a limited set of features. 

ht reference [8]^ the nomialized variance of the magnitude of tiie input baseband 
signal is computed and compared to a predetermined threshold in order to decide in favonr 
of one of the following modulation types: FSK, FM or QAM. Tins method is limited by 
the ijumber of modulation schemes it can identify. As well, it uses only a single feature to 
perform the identification. 

In reference [5], histograms based on the power spectrum of the input signals are 
computed. Freqaencj^ locations and amplitudes are recorded, as well as the locatioii of the 
centre frequency of the overall spectrum. The particular pattenis of each histogram are 
compared to those of tj'pical modulation schemes, such as AM, FSKj PSK or SSB, The 
main drawback of this method is that it uses only the spectral representation of the sigrial 
to perform its computation. 

In refer^ce [10]^ a method is used to discriminate beween an FM aignal and a 
7T/4-DQPSK signal in analog AMPS and digital DAMPS systems. The variation in 
amplitudes of the tu^o different modalation types detenmoes which one is present. The 
main problem v^ith this method is that it is limited to only two modulation types. 

In refei^ence [1 1], the method uses a neural network to demodidate the signal of a 
particular digital cotTsmunication standard. This method differs from the present invention 
in that it identiiies the information content of the signal instead of its format. 

In reference [12], a method is us^d to discriminate between the VSB and QAM 
sigiials that are encount^ed in High Deilnition TV. The main deficiency with this method 
is tliat it is limited to the two modulatiou schemes it can identift'. 



In reference [13], tlie spectral energj'' distribution of the input sdgnal is compared 
with tlic pre-stored energy distributions of FDM/FM signals coBtaiuing specific 
paraEnetora. Recognitioii of a specific foita of signal is declared if the input distribution 
matches one of tbe pre-stored versions. This method is limited to a single form of signal 
feature md to a very specific modulation format, 

III reference [14], the demodulated signal of an FM receiver is olassifled according 
to the voice coding algorithm ttiat processed it This method differs ftom the present 
inv&ntion because classiiicaiioD is applied on the demodulated signal- 
In reference [15], a method is used to generate a decisidfi-tree classifier from a set 
of records. It differs Scorn the present invention in that it does not consider the speeific 
classiflcatioii of modulation formats. 

Rjsference [16] descri-bes a method and ^paratus for detecting and classifyiag 
signals that are the additive combinatioii of a few constant-amplitude sinusoidal 
coDiponents. The main drawback of this schame is that cannot be applied to the 
moduiatioos treated undfer the present invention, except for CW- 

Itz reference [17]^ a sequence of estimated magmtudes is generated from the 
received signal at tJie symbol rate, and the result is compared to 2l predetfennined 
representaticu of Ioqowu voicebpnd digital data rnodem signals. This method is liniitsd to a 
single decision level, as opposed to the series of binary decisions performed under the 
present invention. 

Prior art References 

[1] C.Y. Huang and A. Polydotos, "LilceHhood Methods fer MPSK Modulation 

Classification^ IEEE Transactions on Comniunicatlcns> V ol. 43, Nos. 2, 3 and 4, 
(19S5), 1493 1504. 

[2] E.E. Aszouz and A-K. NandL ^'Automatic Modulation Recognition of 
Communjcations Signals*'^ Klewer Academic Press^ ( 1996), 217. 
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[3] A-K. Nandi E.E. Azzouz^ "Algorithms for AutOTnatic Modulation Recognition 
of CommuiiicatLOTi Signals", IEEE Trainsactioti On Cominunicationg, V ol. 46, No. 
4, (April 1998), 431-436. 

[4] E,E. Azzouz and A^K-Nandi, "Automatic Modulation Recogoition. - 1 & II", J> 
Frankliil Ingiitute, V ol 334B, No 2, (1997), 241-305. 

[5] Robert L. Carrick, William T. Manimig and Robert E. Grimes, US Patent No^ 
4,227,255, October 7, 1980, "Signal Classifier". 

[6] PMlip E- D. Wakeman, US PatetntNo. 4,501,020, Febmaiy 19, 1985, "Spectrum 
Surveillance Receiver System", 

[7] Bruno Lobert and Bruno Scurdillat, US Patent No. 5,27 U036, December 14, 1993, 
*'Methocl and Device f^r the Recognition of Modiilations", 

[8] Nevio Benvenuto, US Patent No. 4,815,137, Marph 21 , 1989, "Voiceband Signal 
deification" - 

[9] Jolin Apostolos and Robert P, Boland, US Patent No. 45166,980, September 4, 
1979:* "Method and Apparatus for Signal Recognition". 

[10] Christoplier Koszarsky, John Nortbcutt and Michael Nowafc, US Patent No . 

5,912,922, June 5, 1999, "Mettiod and Apparatus for Modnlsition Differentiation", 

[11] Alain Chiodim, US Patent No. 5 ,909,675, Jime 1 , 1999, "Device far Recognizing 
Information Conveyed by a Received Signal". 

[12] Ciirl G. Scarpa, US Patent No- 5,636,250, Jime 3, 1 997, "Automatic VSB/QAM 
Modalaticn Recognition Method and Apparatvis". 
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[15] Ronald L. Isaacson and Amy L. Moore-McKee, US PaxeiEit No. 4,845,707, July 4, 

1989, "Frequency Dmsion Mulcipl^x^M Modulation Rjecognition System". 

[14] Chin-Pan Wong and kichard S- YoTOg, US Patent No, 5,651,030, July 22, 1997, 
"Receiver with Signal Classifier". 

[15] R^esli Agrawai and Manish Mdita, Joiuz C. Shafer, US Patent No. 5,870,735, 

Fsbruai^' 9, 1999, "Method and System for Generating a Decision-Tree CiassijRea* in 
Pajrallel in a Multi-Processor System''. 

[16] Neil B. Cox md Edwin L. Froese, US Patent No. 5,353,346, October 4, 1994, 
^^Multi-Frequency Signal Detector and Classifier". 

[17] Nevio Benvenuto and Thomas W, Goeddel, US Patent No. 4,979,21 1, December 18, 

1 990, "Classifier for High Speed Voiceband Digital Modem Signals". 

SUMMARY OF THE INVENTION 

It is now an object of the invention to provide a more straiglitforward and 
computationally simpler meliiod of modulation recognition than neural network based 
methods- 

It is anothier object of tlie invention, to provide a modulation recogijilion method 
whieli directly exploits the fundamental principles of the decision tree based methods in 
order to automatically and accurately perform recognition of a wide varJ.ety of modulation 
formats that are embedded in unkno^-n commuuications signals. These modulation 
formats comprise the foUowing set; CW, AM, FM, FSK, DSB-SC, BPSK, QPSK, -jtJA- 
QP&K, MPSK, NOISE, and OTHERS, 

It is a further object of the present invention to provide a modulation recognition 
method that successfully classifies the signal's embedded modulation amid SNR5 as low 
as 5 dB, earner frequency errors up to +A 50% of the sampling xat^j and carrier phase 
errors Up to 4-/- ISO degrees — all without a j?r/orz knowledge of the syrabol transition 
timing of the signal. 



It is a farther object of the imrontion to extract flic digital complex baseband of 
uijknowxi signals that are measured off-air in real-time, on-line in real-time^ or from 
storage, and then accxjrately classify the signars embedded modul^iion through an orderJy 
series of deoisioE functions- 
It is still a ftixtlier object of the invention to pro-vdde a system for vecogmzmg the 
t^'pe of modulatioii of a modulated signal having a carrier firequency cpmprisixig; 

(a) a receiver section for exiractixig firom the modulated signal a complex baseband 
signal; 

(b) a. pre-classjfication stage for gcaaerating a pre-processod signal fi:om the 
baseband signal; 

(c) means for examimng amplitude variations in the pre-processed signal to 
identify the type of envelope thereof as being one of a constant envelope and 
one of an irregular envelope; 

(d) means for estimating the carrier firequency eicror in the pre-processed sigj^ial; 

(e) means for correcting the pre-processed signal for the carrier frequency errors to 
generate a carrier-corrected signal; 

(f) mean$ for classifying fiie type of modulation in ti^e canier-coTrecled signal^ 
based on the type of envelope identified in the examining step, 

ig) means for classifying the type of modulation embedded in a constant envelope 
signal according to the phase and frequency contents of the pre-processed and 
carrier-cotrected signals^, such modulations being categorized within the set 
{CW, FSK and FM}; and 

(h) means for classifying the type of modnlation embedded in an irregular 

envelope signal according to the phase and amplitude contents of the carrier- 
corrected signal, such modulations being categorized within the set {AM, 
DSB-SC, BPSK, QPSK, ti/A-QPSK, MPSK, OlHER}. 

The method of the invention provides a uniqne decision tree junhitecture that 
automatioalty performs recognition of a wide variety of embedded modulation formats ui 
unkno\vn communicadons signals. The principal method of the invention extracts the 
digital complex baseband of the unknown signal and then determines and classifies. 
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through an orderly series of signal processing functions, the signal's embcdd^ 
modulation scheme with a high degree of accuracy. Without knowledge of the symhol 
transition liming, the method of the inventicfn can be nsed for suocessfiilly perfoimiog 
modulation recognition of measured signals amid SNR's as low as 5 dB, carrier frequency 
errors of up to 4-A50% of the sampling rate^ and carrier phase mora of up to SO 
degrees. The tolerance to carrier frequency errors is limited only hy the frieqnency 
bandwidth of the unknown signal, since the only requirement is to prevent the fireqnency 
eiror from shifting the input signal outside gf the sampling bandwidth of + half of the 
sampling rate. The prefeired method according to die invention is computationally simple, 
making it possible to obsenre the signal during a time frame of preferably less than 100 
msecs- The modalsr arcliitectisne of the mediod of the iuventiou allows for expansion of 
new modulation classifiers. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Exemplary embodiments of the invention will now b« described in more detail 
with reference to the appended figures, in which the referenced munerals designate similar 
parts throughout the figures thereof, and wherein: 

Figure 1 is fi top-level block diagram that illustrates an embodiment of the present 
mvention for measyring and classifying aa off-air unknown signal. 

Figure 2 illnstmtes, in a top-level flow chait, the modulation classification section 

shown in Figure 1 , 

FigTjre 3 illustrates a sample of the spectral power density of the signal and 
embedded noise that can be measured in the pre-cJaasification stage shown m Figure 1 and 
detailed in Figure 2. 

Figure 4 illnstrates ia a graph the degradaiion m the classification success rate of 
four specific signal types versus the normalized cfflicr frequency errors-, without any 
provision for correcting the frequency errors in accordance with this invention. 

Figure 5 is a block diagram illustrating the signal processing fimctions in the error 
correction block 26 of Figure 2 for cocrectrng the residual carder firequency error inherent 
in 0. coxLstant eavelope signal. 

Figure 6 is a flowchart illustratmg liie signal processing fimctions m the classifiea: 
block 30 of Figure 2 for classifying the modulatim of constant envelope signals. 
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Figure 7 is a block diagram illustrating the signal processing fimotions for 
generating the instantaneous frequency components of a constant envelope signal in the 
phase processor "block 42 qf Figtire 6. 

Figure S is a block diagram illustrating tlie signal processing functions in the error 
correction block 27 of Figure 2 for correcting tlie residual carrier firequency exior inherent 
in some irregular envelope signals. 

Figure ? is a graph illustrating the performatics of the signal processing used in 
Figures 5 and 8 for coxrecting the residual carrier frequency error inherent in some 

selected signal? - 

Figure 10 is a flowchart iUusfa-ating the signal processing functions used for 
classifying the modulation of irregular envelope signals in the classifier block 31 of 

Figure 2. i 

Figure 11 is a summary flowchart depictmg the entire modulation classificatian 

I 

process of the present invention, which entails the combination of figures 2, <5 ^nd 10. 

I 

DETAILED DESCmPTION OF THE INVENTION 

Glossary of Acronyms | 

For convenience, a glossary of acronyms used m the description of ttie present 

invention is given bfelow: 



AM 


Amplitude Modulation 


AMPS 


Advanced Mobile Phone Services 


ASK 


AmpHtnde Shift Keying 


AWGTvI 


Additive "White Gaussian Noise 


BPSK 


Binary Phase Shift Keying 


BT 


Product of a filter band-width B and the symbol period T 


BW 


Bandwidth 


CW 


Continuous Wave 


DSB 


Double Sideband 


DSB-SC 


Double Sideband Suppressed Carrier 


FFT 


Fast Fouxieor Transform 



FSK Frequency Shift Keying 

FM Frequency Modulation 

GMSK Gaussian Miniaium Shift KeyTUg 

GSM Global System for Mobile (Pan Europeau Digital Cellular) 

IS-54 Interim Standard #54 of the Tfilecommimications Txidustry 

Assodation (also known as US Digital CelMar) 

MPSK M-ary PSK 

MUX Multi-Clianael or Multiplexing 

PSD Power Spectol Density 

PSK Phase Shift Keying 

QAM Quadrature AcapUtude Modulation 

QPSK Quattemaiy Phase Shift Keying 

71/4-QPSK 7i/4-shifled Quaternary Phs^e Shift Keying 

RMS Root Mean Squ^ 

SNR Signal to Noise Ratio 

SSB Single Side Band 

VCO Voltaga Control OsciUator 

VSB Vesti^al Sideband 



Termindlogy 

Throughout tlie description of the present invention the follomng terms are: used in 
accordance with their reapeetLve definitions given below: 

Linear vector' quantization '""Vector quantization" is a proces3 by whicli a "vectorj, 

defined over a continuous multi-droiensional vector space, is 
assigned one of a limited number of pro-determined 
r^resenlative vwtors. Quantization occuors at the vector 
level because &d infinite number of vector possibilities is 
mapped onto an Msemble containing a finite number of 
vectors. The mapping is peifbtmedby dividing the originaJ 
vector space into a number of regions, each of which 
corresponding to one of the finite number of pre-determined 
representative vectors. "Learning vector quantization" is the 
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process that perfoims th^ adaptive segmentation of the 
infinite vector space in the pre-determined regiODS dormg a 
traioirig period or during normal operaliaa of tlie quantizer. 



Training set The trainitig set is an ensemble of known vectois, signals, or 

other inputs provided to a system that estimates ceit&in 
parameters over the inputs. The system^s internal state is 
thus allowed to adapt to the training set so that it can 
thereafter perform parametric estimations on tanknown 

inputs. 

C^TitTQid frequency The oentroid freqxiency of a given pgwer spectral deasity is 

the frequency correspondijig tg the "cwtre of ma^s" of such 
power 5pectr?Ll density. It is a purely geomelxic propearty of a 
given density i$ expressed as: 




where the d-s axe the Discrete Fourier Transform (DFT) 
coefficients of the input signal and die are the 
corresponding fiequencies. 



Zero padding In the computation of the Discrete Fourier Transform (DFT) 

of a given discrete-time signal vector x(n) of finite length 
zero padding refers to the appending (the padding) of a 
certain number L of samples (all equal to :?:em) at the end of 
tbe vqctpr x(n)- The DFT then computed on a new 

discrete-titne signal vector [x(n) 0 0 0 of length. NH-L. 
This process has the effect of increasing the frequency 
resolution of the origioai TSf-poiat DFT, since it is then 
computed for N-4-L different frequencies over the same 
frequency range. 
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Unwmpped phase Unwrapping the phase of a complex signal vectQi refers to 

the process of making the corresponding Ttiodulo-2K phase 
vector continuous over the range from minus inflnity to plus 
infinity, by adding niiilliples of + 0T-27r as needed. 

Digital complex baseband Tho complex baseband representation of signals (also 

refcrced to as "complex envelope rcjjresentation") serves the 
purpose of mathematically expressing a real bandpass signal 
in a format which contains the amplitude and phase 
information of the signal but not its carrier or reference 
freqviency. This representition is an extension of the familiar 
two-dameixsional phasor (vector) representation of sinusoidal 
signals. 
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Consider a real bandpass signal with, a narrow bandwidth, 
concentrated 'iround a. reference or carrier froqumcyfc^ This 
real signal can be matliematjcally expressed as: 

where a(t) is the signal amplitude and <p(t) is the tim&- 
vaiying phase. Tlie modulation of the signal is contained in 
its amplitude and its phase (the fireqaency modulation is 
equal to the time derivative of the phase). The signals a(£) 
and ^t) show all the -modulation mfonnation earned by sft), 
and must be preserved in the complex baseband 
representation. The ^bove equation cm also be expressed as: 

Substituting: 

%ve then have 3. signal: 

u(t) = jc(/) + JyiO 

which is the complex bas^and equivalent of s{t) and 
contains aU the arr^litude and phase information of the real 
bandpass signal. 
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Overvie-w of Pr^ent Inrentioii 

Figure 1 is a top-level blocQc diagmm that shows one aspect of the present 
inventtoiT - the measurement of off-adr unknown signals. This figure includes an antenna 
1, a receiver section 2 and a modulation classification section 3, The receiver section 
inclades an amplifier and filter 4, a ftequency down-converter 5, a. signal sampling 
processor 6, and an analog to digital (A/D) converter 7. Each of these receiver circiiits are 
generally recognized in the prior ait. The xmknown signal, which may comprise a plurality 
of channels, is collected ofF-air at the antenna 1 and then amplified and filtered hy step 4. 
This function also server to i$olate the channel containing the desired cairier ficquencj'- 
(/g). The filtered signal is then doTi^TL-converted in firequmcy in step 5 to the analog version 
of the complex ba^ebfmd signal. While this process rednces /g to near zero Hz, residnal 
carrier frequency errors still remain within the converted signal and must be dealt with 
prior to the classification process in order to ensure that the modulation format is 
accurately recognized. From the do^Ti-converter 5, die analog complex baseband signal is 
sampled by step 6 at regular time intervals of T^ = l/F^ seconds, and the resuLting time 
samples are quantized by the A/D converter 7 to produce a $eries of quantized time 
Simples known as the digital complex baseband signal Quantization is performed within a 
finite number of levels corresponding to the number of bits of the A/D converter- Hie 
digital complex baseband signal 8 is then fcrwarded to the modulation classification 
section 3y the principal embodiment of this invention, where the digital signal processing 
for classifying the modulation schemes begins. 

The modulation classifier 3 of Figure 1 illustrates a simplified view of the 
modulation classification process, consisting of a pre-classification sta^e 9, a signal 
envelope decision test lO^ and two separate classifLcaiion stages - one for processing 
signals with a constant envelope (11 and 13} and the other for processing signals with an 
irregular envelope ( 12 and 14). The purpose of the pre-classificatiou sta^e 9 is to (a) 
determine if sufficient signal pdwer is present, (b) estimate the signal's bandwidth and 
Qentroid frequency^ (c) perform the first level of cairier frequency error correcliorn^ and (d) 
filter the out-of-band background noise power. The first classification test is performed in 
decision step 10 where constant envelope signals ai^e discriminated jSrom those of irregular 
shapes. If the signal has a constant envelope, its modulation format is categorized v^ifhin 
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the set {CW, FSK, FM}, Conversely, if significant envelope variations are detected, the 
signal's modulation format is categorized wittrni the set {AM^ DSB-SC, BPSK, QPSK, 
7C/4-QPSK, MPSKj OTHER}. Followixig this test, the residual carrier frequency errors 
superimposed on either signal are cancelled when required, in the eiror-correction steps 11 
and 12 respectively, and the processes for examining, recognising and classifying the 
modnlattoii ficheraes on the respective constant envelops or irregular envelope signals are 
perfoiMed in the classifiers 13 and 14, The phase and frequency content of the constant 
envelope signal is examined in classifier 13 in order to classify one of the {CW, FSK, 
FM} modulation formats. The phase and amplitude content of the iTregul^r envelope 
signal is examined in classifier 14 in order to classify one of the {AM, DSB-SC, EPSE, 
QPSK, 7t/4-QPSK, MPSK, OTHER} modulation formats. 

Pre-classification Stage 

Figure 2 shows m a flow chart a more detailed view of the Modulation 
Classificatioti Section 3. The process begins at input IS where a block of digital complex 
ga^xiples is obtained firom the output S of the Receiver Section of Figure 1. These samples 
3re forwarded to a signal detector 16 where Fast Fourier processdng determines if 
sufficient signal power is available to establish the presence of a signal. This process is 
performed with the aid of an external noise power module 17 that establishes a power 
threshold against which the power spectral density of the observed signal is oompaied- 
The noise power 17 is also used in by tiie BW and centroid estimation process 20 to set 
another toeshold for computing fhc centroid frequency and estimating the bandwidth of 
the signaL As iUu$trated in Figure 3 j the noi&e power 17 corresponds to the level Pv^ . A 
signal is djeclared present at decision step IS if it bears some Jxeqtiency components above 
a threshold (P^ h- Z) dB, where Z is a given value relating to the probability that a noise 
jSrequency component has a power level above the (P^ H- Z) dB threshold. Thus, if the 
signal energy is not above the noise-relaied thresihold, the SNR is coiusidered too small to 
allow further pniicessing and the signal is classified as NOISE 19. However^ if a signal is 
declared present, then the gross estimate of the carrier frequency inherit in the signal is 
determined in die B W an.d centroid estimation ^ep 20. This operation is done by 
computing the centroid of the portion of the power specti'al density that is above a certain 
threshold (Pw + Y) dB, as shown in Fi,gure 3. This threshold is lower than the one used in 
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Step 18 to detect the signal presence. Using Fast Fourier processmg, the centroid frequency 
is given by: 

f=r^ (1> 



5 where PSD[i\ is tiie power spectral daisity value coirespoiidiDg to frequaicy^, {md the 
summation is computed for the frequencies corresponding to power spectral density values 
PSD[f] sibQve me threshold (P^, + Y) dB. Frequency translation of the signal is then 
peiformed in the gross error correction process 21 where the carrier frequency is coirected 
by an amount equal to the centroid drequenc>'. The frequency bandwidth is also estimated 

10 m step 20. This operation is done by Heiecting the bandwidth corresponding to a pre- 
selected percentage of the sum of the power spectral dcnsit>' values PSD[i] above the 
given threshold (P,- + Y) dB. This estijiiated bandwidth is then used in step 22 to filter the 
out-of-band noise power firom the fi^quency-translated signal in step 20. The fitter is 
selected fhini a bank of pre-stored filters. 

1 5 The process for detecting aud classifying the modulation format embedded in the 

signal is performed in $tep3 23^, 30 md 31. Rapid and precise classification is possible 
provided that the residual caixier ftequency acror mherent on the output sigads 28 and 29 
is kept to witiiin 0.001% of the sampling freqamcy Fs used hi sampling step 6 of Figure L 
This limitation is illustrated in Figure 4, where the success rate for four different 

20 modulation types is maintained at the zero-offset levels for frequency errors up to 10'=^ x 
Ffi. It should be noted, however, that the presence of amplitude variation on the complex 
baseband signal is insensitive to any phase vector rotation caused by the inherent carrier 
fi-equency error. The present invention takes advantage of this feature when discximinating 
between constant envelope signals and irregular envelope signals in decision test 23. This 

25 binary classification test is performed by using the maximum o f the squared Fourier 
Transform on the normalized signal amplitude, given by 



]DFr(a)]" . 
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where /is the j&eq.uency, DFT(.) is tte Discrete Fourier TrarLsfomi (computed fiom tlie 
Fast Fourier Transform), Ns is the number of samples in the input block, and a is the 
amplitude vector centred on zero and normalrzed by its raean. Mathematically, vector a is 
expressed as 

'TTI 

a = 1 

where x is the observed filtered vector^ 5c is the vector of magrdtudjes for the elements of 
vector Xj B[\x\] is the average of the magnitude elements aver the vector il. The test 

paiaineter /max is totally mdependent of the phajse variation in x and is therefore 
unaffected by residual carrier Jrequency errors. At the same time, it allows for the 
identification of information in the signal envelope, which^ in most cases of interest, 
results in the rehab J e discrimination between constant envelope tod irregular envelope 
signals. The peiformancs obtained in Equation 2 is very good for SNRs as low as 5 dB, 
where the SNR is defined over the sampling bandwidth. 

Frequency Error Correction of Constant Amplitude Signals 

When test 23 determines that the signal has a constant envelope, the process moves 
to the error-correction step 26 where the residual carrier frequency error is estimated and 
corrected. The puipose of step 26 is to caooel out the carrier frequeoncy earor inherent in 
the constant amplitude signal 24 an order to detennine if the signal is a continuous wave 
(CW or pure tone). Recognition and classification of a QW signal cannot be performed on 
the constant envelope signal unless its inherent carrier frequency error has been corrected. 
The operation is performed in the frequeDcy domain as shown in Figure 5, consisting of 
three signal processing operations: a coarse estimation of the earner frequedocy ( 32, 33, 34 
and 35); a refined estimation of the earner frequency (37); and the carrier frequency error 
correction, of the constant envelope signal ( 39 and 40). The process begins at st^ 32 
where the normalized squared signal is computed. The result is then processed by a Fast 
Fourier Transform (FFT) 33 with zero padding^, squared in step 34, and forwanied to a 
coarse frequency estimatoo: 35 where a search is done to locate the frequency line with 
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majiimiiin powet. Note that both the zero padding and the squaring of the output of the 
FFT improve the frequency resolution of Oie coarse frequency estimator 35, The output 
36^ denoted hyj\um is forsvarded to a fixie firequeiicy estimatOT 37 wliere the reEned 
estimation of the carrier frequency is performed. This signal processor (37) 5l1so receives 

the constant envelope signal 24 and computes the estimated residual carrier frequency 

^Ttcr A/" at 38, as follows: 



where PSDDj/] is the derivative of tho power spectral densitj"" at frequency^ /] is a 
frequency given by 



where Fs is th^ sampling jBrequency and the number of samples used in the FFT (that is, 
including zera padding). The plus sign is s&Iccted for fi if PSDDj)^^] is positiv^^ and die 
negative §ign is selected if PSDDj)^] is negative, .PSDD[,] is obtained as 



where Re[_] indicates the real part, S(f) is the Discr^e Fourier Transform (DFT) of the 
filtered observed signal at 24, S'^^fi is its complex conjugate, and S'{f) is its first derivative 
^ntb respect to frequency, given by 



where xjc is a sample of the filtered observed signal vector at 24 and N is the nmnber of 
samples in the input vector, 

The ^ttttiated firequency error Af at 38 is used to set the firequency of a digital 
VCO 39, The output of tM? VCO is then multiplied by the constant envelope signal 24 in 
a multiplier 40 to cancel out the signal's residual carri^ frequency error. 



A/ 



(3) 



PSDDt/^]-PSDD[/J 



PSDD[/] = 2Re[5^X/)5'^(/)] 
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Classification of Constant Envelope Sign^;? 

Tlie output 26 of Figure 5 (and Fi^e 2) is a frequency-translated, cotistant 
envelope signal, rhis sigoal is now forwarded to classifier 30 of Figure 2 where its 
embedded modalation is fbrther examined to deteimine if it is CW, FM or FSK, This 
classificatioTi process is illustrated in Figiire 6 and consists of two steps: (a) a test for 
signiacant plmse variation to detemnne if fhe signal is CW or ftequency modulated ( 41), 
and (b) a series of tests for determining specific- frequency ctaractexistics to ena'ble 
discrimination between FM and FSK formats (42, 43, 44, and 46). Tho process begins at 
decision step 41 where the frequency-translated, constant envelope signal 28 is examined 
to deteTiniiie if significant phase variations are present vA&m its envelope. Since CW 
signals have a near-zero variance in the unwrapped phase (direct phase) and frequency- 
modulated signals have relatively high variances, a test that deteimmes the variability of 
the unwrapped phase is used. However, tliis unwrapped phase test requires care when fhe 
signal vector amplitude is small because phase variations increase as tiie SNR decreases. 
The difficulty of ihis problem is mitigated by discarding any data sample whose amplitude 
is below a threshold. Such threshold lowers ths probability of residual pliase variance in 
the CW signal, while not affecting the phase variance of frequency-modulated sisals. An 
amplitnde threshold equal to tli€ meao of the amplitudes of the sampled signal is chosen 
under the preferred embodiment of this invention. Di5criminatioa between CW and 
frequency modnlated signals is obtained by first retaining the samples of the signal 2S at 
an amplitude above the amplitude threshold, and then comparing the value of the 
unwrapped phase of these stoiples to a phase threshold whose value is generally optimized 
through pgjnpnter simnlation experiments. If little or no phase variation exists, decision 
test 41 determines that the signal is CW and notifies the phase processor 42 not to proceed 
with further classifloation. Conversely, if significant phase variations occur in fhe signal, 
step 41 declares Ihe presence of fiequ^cy modulation and notifies the phase processor 42 
to proceed with the next stage of classification process for discrimination between FM or 
FSKfoimats. 

To perform the classification between FM and FSK formats, the frequency- 
nncorrected, constant envelope signal 24 is used instead of fhe frequency-translated signal 



19 



28 because the process for oompating tJie latter signal (described in Figure 5) may create 
large frequency errors for non-CW signals, and such large errors would create problems in 
detenniuing the FM and FSK classifications. On the other hand, the residual carrier 
frequency error inherent in signal 24 is sufficiently small to allow for accurate 
discrimination between FM and FSK signals. In order to perform such discrimination 
be1>yeen ¥M and FSK, the kurtosis coefficient of the instantaneous frequency content of 
the signal is need. This coefficient represents a measure of the flatness of the signal 
frequency distribution, and its value tbr analog FM is usually different from that of digital 
frequency modulated sisnals- The kirrfcosis coefScient also has a goodbehaxdour at low 
SNRs. The coefficient is given by: 



where^tt) is the instantaneous frequency (about the mean) at time and E[/?y] is the time 
average of/ft) computed over the length of the observed si^^ml. 

The process begins at the phase processor 42 where the instantaneous frequency is 
obtained by computing the phase derivative of the constant envelope signal 24. This 
computation is further illustrated in Figure 7, consi^g of signal processing steps 49, 50, 
51, 52, and 53. The phase of the constant envelope signal 24 is computed in step 49 to 
produce the instaaitaneous phase signal 50, The instantaneous phase si^al is then 
processed througli a low-pass filter (51 and 52) in order to hmit the noise enhancement 
effects of tlie derivative fimction in step 53. A simple computation in BW estiEnatioo 51 
roughly estimates the edSective bandwidth of the phase sigp.al which, an tura^ selects fifim a 
library one of a small number of low-pass filters that has a cut-off frequency slightly 
higher than the signal bandwidth. The output of the dmvative function 53 is the estimate 
of the instantaneous frequency 54. This estimated value is required in the computation of 
the kurtosis coefficient in step 44 of Figure 6, It is also required to detect the possible 
presence of a low-modulated commercial FM signal in decision step 43. The FM 
modulation of such a signal contains the presence of silent periods, and, when 
demodulated, it produces a single discrete tone. Other forms of FM signals also contam 
such tones, such as the analog AMPS signals. Thus, the test perfonned on the constant 
envelope signal 24 in step 43 serves the purpose of detecting the presence of such a tone. 
If the tone is present, the signal is declared to be FM (commercial or AMPS)- Conversely, 
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if such tone is not present, the signal is ciHier narrowband FM or FSK. Following this 
process^ the kurtosis coefficient is tested. FM signals usualty praduce alairtosi$ coefficient 
that is larger than 2.5. Converselyj the coefficient prodncsd by FSK signals i? smaller 
because such signals are mqre compactly distributed. Thus, if the kurtosis coefficient is 
below the 2,5 threshold, a decision is made in favour of FSK. When the loirtogis 
coefficient is higher than the 2,5 threshold, an mitial decision is made at 45 in favour of 
FM. However, a further classification test is required at this point because FSK signals ' 
v/ixh highly filtered data (sMch as GMSK signals) may also have a high kurtosis 
coefScient- A frequency-modulated signal whose digital modulation contains a symbol 
rate off^y^ and a frequency deviation of a multiple otf^JA^ produces discrete firequency 
lines in the power spectral density (FFT) of the squared signal. For examplSj MSK and 
GMSK signals have a frequency deviation of 2 x tymM -fsyn^ /2 and a modulation level 
M = 2. The FFT of the square of either signal produces two discrete frequency lines. 
Accordingly^ the discrimfnatjon process is further refined in decision test 46 where the 
FFT of the squared constant envelope signal 24 is computed. (This computation is the 
same the one done in steps 33 and 34 of Figure 5). Thus, if two or more discrete 
frequency lines are delected, the signal is classified as FSK at 48* Conversely^ if decision 
test 46 produces less than two lines^ the signal is classified as FM. 

This concludes the detailed descriprion for the classification tests performed on 
constant envelope signals, 

Frequeacj^ Error Carrection of Irregular Envelope Signals 

The following stq5S describe the classification processes for signals that have 
irregular envelopes. When decision step 23 of Figure 2 detemiines that the measured 
signal has au iixegular envelope, the process moves to the eiTor-correction processor 27 
where the residual carrier frequency enpors are estimated andj corrected. This operation is 
performed in the frequency domain as illustrated iji Figure S] and is similar to the 
operation described in Figure 5. Again, the process consists of tfciree signal processing 
steps: a coarse estimation of the carrier firequency (55, 57, 59^ and 6t); a refined 
estimation of the carrier frequ^cy (62]; and the carrier frequency error correction of the 
irregular envelope signal 25 (65 and 66). The process begins at step 55 where the square of 
the irregular envelope signal 25 is performed. The resulting signal is then processed by 
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the FFT 57 with zero paddiag^ squared in step 59, and foiwarded to a coarse frequency 
estimator 61 where a search is done to locate the freqiieEiicy line with maximum power. 
Zero padding and the squaring of the EFT output 58 perform the same resolution ihnctions 
as those described in Figure 5. Depending on the modulation format of the iuilcnov»'n 
signal, a peak: in tJi^ firequeiiGy will appear in path 50- Choosing the frequency Hne with 
the highest peak, the process moves to a fine frequency estimator 62^ where the carrier 
frequejicy estimation is ftuther refined^ This signal processor (62) also receives the 
norniafeed squared version 56 of the irregular euvelorpe signal 25 and, using the one-step 
of the secant method as described in equation 3, computes the estimated carrier frequency 
error. This value is halved In step 63 according to the squaring non-hnearity ih step 55, 
The result Af at 64 is used to set the frequency of a digital VCO 65 whose output {exp(- 
j27r4/t)} is then multiplied by the irregular envelope signal 25 in a multiplier 66 to cancel 
out the signal's residual carrier freqaeaey error. As illustrated in Figure 9, the foregoing 
frequency error coixectioa process is enable of obtaiuing a normalized RMS frequency 
estimation error lower than 10"^ for the SNRs of interest (ie., larger than 5 dB). This carrier 
recovery method is therefore sufficiently precise to avoid performance degradation due to 
the limitations described in Figure 4. However, it should fae noted that the frequency error 
estimation scheme performed in Figure 8 does not produce ati adequate estimation of the 
carrier firequency for MPSK signaJs. This limitation notwithstandiog^ the follow-on 
process leading to the classification of MPSK uses the irregular envelope signal with 
residual carrier frequency errors at 25 of Figure 2. 

Classification of Irregular Envelope Signals 

The frequency-translated, irregular envelope signal 29 of Figure 8 (and Figure 2) is now 
forwarded to classifier 31 of Figure 2 where it is furfheor ttiamincd to detemiiae if it is aiuougst th.e 
set {AM, DSB-SC, BPSK, QFSK, MPSK, OTHER). Ibis classiacation process is illustrated in 
Figure 10 and comprises the following signal processing steps: discriittmatioiL between one- 
dimensional and t^vo-dimensional signals (67); classification of one-dimensional signals (6S and 
70); and classificatioii of two-dmensioTial signals (73, 75, and 77). The process begins at decision 
test 67 where the firequency-translated, ircegular envelope gjgtml 29 is examiued to debenuine if 
any modulation uiformation is carried In the phase of the signal. If such informatioii is not present, 
the signal modxU^tian corresponds to the actual baseband, such as AM (transmitted carrier), DSB- 
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sc. or EPSKj and is recognized as one-dimcnsionaL These signals recognised from their 
absolute ccTitied un^vrapped phase sequence which has a small -variance, practice, the -process 
requires phase unwapping because the initial carrier phase is random. However, pjiase 
unwrapping is Ycry sensitive to notse^ paiticalarly on DSB-SC and B?SK signals due to theix 1 SO"" 
5 phase transitionib. Moreover, aB.y phase unwrapping error will seriously Lacrease the variance of the 
absolute centred phztse. To avoid phase unwrapping of DSB-SC and BPSK signals, a new quantity 
is introduced in step 67 Galled the "absolute" phase, and is defined as 

tio = z(\m\+ j\Q(t)\) 

10 where I{f) and Q{t) are the mphase and quadrature samples at tiitie ^ By tekiag £he 

absplu-te values of the real and iniaginaiy parts of a signal, a phase sequence between 0 
and 7^'2 radians is produced. The phase sequence allows for a significajiit reduction in the 
size of the observation space nequixed to make a decision. For DSB-SC signal:^ the 
variance of the phase of the resulting signal is generally smalL Sniall variance is also taie 

15 for BPSK and AM signals. For two-dimensional signals beariTig some phase infortnation 
(such as QPSK mdMFSK signals), the variance of the absolute phase is high^ tending 
towards the variance of i uniformly distributed random variable in die interval [0, 
radians. To pro\ddje a better separation between one and two -dimejisional baseband 
signals, a threshold oii the amplitude of the signal is nsed in step 67, and the samples 

20 below this threshold are discarded. The use of the amplitade threshold reduces the 
variance of the phase on one-dimeiisional signals, Tvithont significantly affecting the 
variance on t^^'^o-dimensional signals. Simulations show that the best results occur when 
the said threshold is set equal to the mean of the ampHtide of the phase sequence. If 
decision stiep 67 determines that the signal modulation is one-ditnensionaL the process 

25 nioves to decision test 63 wjiere the variTOce of the unwrapped ph^e (direct phase) is 
conaputed. Due to the presence of n radian jumps in the instantaneous phase of DSB-SC 
and BPSK signals, the said variance is much higher for these signals thim for AM signals. 
Phase unwrapping is useful in this case because errors of this quantity are unlikely for AM 
signals, thus providing a good discriminating feature. In order to reduce the phase vatriance 

30 for AM signals, a threshold equal to the mean amplitude is set on the amplitude of die 

signal samples in step 68. An indication of low phase variance thus classifies the signal 9.5 
AM at 69. Note also that ASK modulation is a digital form of amplitude modulation. ASK 
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signals are therefore classified as AM signals at 69. If the pha5o variance is highj decision 
step 68 determines that the signal is not AM and moves to decision step 70 where a fuxtiier 
test H made to deteimine if the signal is BPSK or DSB-SC, Unlike BPSK signals, DSB- 
SC signals bear an envelope for ^^^hich. the amplitude varies substantially overtime. 
5 Accordingly, step 70 computes the variance of the envelope and uses it to discrimttiate 
between a DSE-SC signal at 72 and a BPSK sigaEil at 71, 

Retaining to step 67» if the test determines that the signal is two-dimensional, the 
process moves to decision steps 73, 75 and 77 where PSK signals are separated from 
QAM and other unidentified modulation types. To perform the classification between 

10 QPSK, MPSK and OTHER types of two-dimensional formats, the firequency- 

imcorrectedj constant envelope signal 35 fiom decision step 23 is used instead of the 
fireqnency-translated signal 29 computed in Figure 8, because the residual caixiex 
frequency error inherent in signal 25 is sufQciently small to allow for accurate 
discrimination between the QPSK, 7I/4-QPSK, MPSK and OTHER signal?. To initiate this 

15 classification in decision step 73^ a simple test on the variance of the signal amplituds is 
performed. If the variance of the iimplitude is below a given threshold^ the signal is 
assumed to be a PSK signal. Even thaugh band-limited PSK signals do exhibit amplitiide 
variations, such variations are much less noticeable than those for QAM, SSB or VSB 
signals. This test is therefore sufficient to discriminate PSK signals from most of the other 

20 two-dimensional signals as OTHER at 74, If iJae test deteimmes that the signal is PSK, the 
process moves to decision step 75 where the signal is fhrther examined to determine if its 
embedded modulation is QPSK, 7iy4-QPSK or MPSK, This test as performed by 
computing the fourth power of the sigaal 25. QPSK signals produce a single peak in the 
power spectral density of the resulting signal, while ttM-QPSK signals produce two peaks 

25 separated in frequency by twice the baud rate. Therefore^^ if the test in step 75 results in 
One peat^ the signal is classified as QPSK at 76. Otherwise the process moves to decision 
test 77 where the combination of two peaks classiiies the modulation as 7c/4-QPSK at 78 
and the absence of any peak classifies the modulation as MPSK at 79. 

This concludes the detailed description &ir the classification tests of irregular 

30 envelope signals. The description covering the overall decision tree process as illustrated 
in Figure 1 1 , is given below. 
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Overall Description of the Preferred EmbodiBieilt 



1 



p-gyj-p 1 1 illustrates a complete overview of tiue modulatLon classificatioii section 3 
of Figure 1 and collects togellier the signal processing fonCtiotls illustrated in Figures 2, 6 
5 and 10. The process begins at input 15 where a block of digital compleK samples is 
obtained jfrom the output S of the Receiver Section of Figure I. These samples are 

r 

forwarded to a signal detector 16 wli(sre Fast Fourier pTccessing determines if suffioiem 
signal power is available to establish tlie presmce of a sigrxaL This process is pcaformed 
with the aid of an external noise power module 1 7 that estabHshes a powa: tiireshold 

IQ against which the power spectral density of the observed signal is compared The noise 
power module 17 is abo used in step 20 to set another threshold for computing the 
centroid frequency and estimating the signal bandwidth of the signal. If the signal energy 
is not above a noise-related threshold, the SNR is cornsidered too small to allow further 
processing and the signal is classified as NOISE in step 19, However, if a signal is 

15 declared present, then the gxoss estimate of the carrier fiequency inherent m the signal is 

determined in step 20. This operation is done by computing liie centroid of the portion of 
tixe powet spectral density that is above a certatn threshold. Frequency translattoi? of the 
signal t$ then performed in correction step 21 where the carrier firequency is corrected by 
an amount equal to the centToid frequency. The frequency bandwidth is also estimated in 

20 estimation step 20. This estimated bandwidtli is then used in filter step 22 to filter the out- 
of-band noise power &oni the ftequency-traixslated signal in estinaadon step J<l. The 
process for deteetitig and classifying the modulation founat embedded in the signal starts 
in decision test 23 where discrinunation between constant envelope signals and irregular 
envelope signals is performed- Wlien decision test 23 determines that the signal has a 

25 constant envelope, the process moves to the error-correction step 26 where the residual 
earner jft-equency error is estimated and corrected. The classification of constant-envelope 
simals starts in decision test 41 and consists of ttie foUowng steps: (a) a test for 
significant phase variation to detennine if the signal is CW or firequency modulated (41), 
aad (h) a series of tests fox detetmining specific &equency diaracteristics to enable 

30 discriTimiation between FM and FSK formats (42, 43, 44 and 46), If little or no phase 
variation exists in the frequency-translated signal, decision test 41 determines that the 
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signal is CW and notifies the phase processor 42 not to proceed with further classij&catioTi, 
Conversely, if significant phase variations occur in the signal, decision, test 41 declares the 
presence of frequency modnlatioii and notifies the phase processor 42 to proceed with the 
next stage of classification process for discrimiBation betw'^een FM or FSK formats. To 
5 perform the classification between FM and FSK foraiats, the frequency-uncorreoted, 
constant-envelope signal 24 fi"oni decision step 23 is tised. The process begins at 
computing step 42 where the instantaneous frequency is obtained by computing the phase 
derivative of the constant envelope signal 24, This estimated value is required in the 
computation of the kurtosis coefficient in decision test 44. It is also required to detect the 

10 possible presence of a low-modulated commercial FM signal in decision test 43. The FM 
modulation of such a signal contains the presence of silent periods, and, wiien 
demodulated^ it produces a single discrete tone. Other forms of FM signals also contain 
such toneSp such as the analog AMPS signals. ThuSj. the test performed on the constant 
envelope signal 24 m decision test 43 serves the purpose of detecting the presence of such 

15 a tone. If the tone is presenij the signal is declared to he FM (commercial or AMPS). 
Conversely, if such tone is not present, the signal is either narrowband PM or FSK. 
Following this process, the kortosis coefficient is tested in decision test 44. FM signals 
usually produce a kurtosis coefficient that is larger than 2.5. Conversely^ the coefficient 
produced fay FSK signals is snoialler because such signals are more compactly distributed- 

20 Thus, if the kurtosis coeflicient is bdow the 2,5 threshold, a decision is made in favour of 
'FSK:. ^Tien the kurtosis coefficient is higher than the 2.5 threshold^ an initial decision is 
made in favour of FM. However, a further classification test is required at this point 
because FSK sig^ials with highly filtered data (such as GMSK signals) may also have a 
high kurtosis coofficient. The discTimiitation process is finther refined in decision step 46 

25 where the FFT of the squared constant envelope signal 24 is computed- If two or more 
discrete frequency lines are detected, the signal is classified as FSK. The estLmafce of the 
n Limber of symbols Mis obtained in coimtlAg step 47 by counting number of FFT 
peaks. Conversely, if decision test 46 produces less than two lines, the signal is classified 
asFM. 

30 WlLcn decision step 23 detemoines that the measured sig;nal has an irre^lgqc 

envelope^ the process moves to the error-correction step 27 where the residual carrier 
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frequency errors are estimated and corrected. The frequency-txanslated, irregular envelope 
si^al output of error correction step 27 is now forwarded to decision test 67 where it is 
fortlier examiaed to determinfe if it is amoagst liie set {AJVL, DSB-SQ BPSK, QPSK, 
MPSK, OTHER}. This classification process comprises of the followiag signal processing 
steps: discrimination between one-dimensional and two-dimensional signals (67); 
classification of oii&-dimensional signals (68 and 70); and. classiEcation of two- 
dimensional signals (73, 75, and 77). The process begins at decision test 67 where the 
frequency-translated, irregular envelope signal is examined to deterffline if aa^'^ modulation 
information is carried in the phase of the signal. If such information is not present, the 
signal modulation corresponds to the actual baseband, such as AM (txansmitted carrier)^ 
DSB-SC, or EPSK, and is recognized as one-dimensional. If decision test 67 determines 
that the signal modulation is one-dimensional^ the process moves to decision test 68 where 
the variaaoe of the unwrapped phase (direct phase) is computed. An indication of low 
pha5? variance thus classifies the signal as AM. If the phase variance is high, decision test 
68 determines that the signal is not AM and moves to decision step 70 where a further test 
is made to determine if the signal is BPSK or DSB-SC- UnUke BP5K si^jials, DSB-SC 
signals bear an envelope for which the amplitude varies substantially over time. 
Accordingly, decision test 70 computes the variance of the envelope and uses it to 
discrimmate between a DSB-SC signal from a BPSK signal. 

Retuming to decision step 67, if the test determines that the signal is two- 
dimensional, the process moves to decision steps 73, 75 and 77 where PSK signals arc 
separated from QAM and other imidentified modulation types. To perform the 
classification between QPSK, 7t/4-QPSK, MPSK and OTHER types of two-dimensional 
formats, the frequency-uncorrectcd, constant envelope signal 25 from decision test 23 is 
nsed. To initiate this classification in decision step 73, a simple test on the variance of the 
signal amplitude is performed. If the variance of the amplitiide is below a given threshold, 

the signal is assumed to be a PSK signal. Otherwise, it is classified as OTHER. If the test 
determines that the signal is PSK, the process moves to decision step 75 where the signal 
is further examined to determine if hs embedded modulation is QPSK^ tiM-QPSK or 
K'IPSEl This test is performed by computing the fourth power of the signal 25. QPSK 
signals prodace a single peak in the power spectral djensity of the resulting signal^ while 
7T/4-QPSK sisals produce two peaks separated in frequency by twice the band rate. 
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Therefore^ if the test in step 75 results in one peak, ftue signal is classilSed as QPSK, 
otliearwise tlie process moves ta decision step 71 where the combination of two peaks 
classifies the ixiodulation as jt/4-QPSK and the absence pf any peak classifies the 
modulation as MPSK- 

5 Different aspects of the embodiment in Figure 1 axe possible. For example^ the 

signal coming to the Receiver Section 2 can be measnred from a tapped wireliiie, or &om 
another receiver or pre- amplification section. The signal may also have been amplified and 
filtered fay an alternate receivsr^ and he introduQed before the firequency down conversion 
5. Similarly, the down conversion to complex base^band eauld be perfoixaed on an 

10 alternate set of hardware^ and the signal could be iatrodaced before the sampling fbnction 
Furthermore, the sampled digital complex baseband signal 8 could be obtained from a 
receiver section with a different configurations c>r it could be retrieved j3-om a digital 
memory location where it had been previously stored. Note that, for practical reasons, the 
sampling function 6 and the A/D converter 7 are usually implemented in a single unit, and 

15 that it is not very likely that the signal could be introduced before the digital conversion 
function 7. Such a procedure is nevertheless coiuceptually also possible. 

PERFOllMANCE EVALUATION 

In order to demotistrats the methods described in the present hiveDtion, 500 
20 simulated signals of each of the modulation types have been generated and processed. 
These signals covered a wide variety of parameters as described in tiie next section, A 
sampling frequency of 48 kHz was used covering a bandwidth slighdy larger than the 
occupied bandwidth of most narrowband connnunications signals. Sequences of S5,3 msec 
(representing 4096 samples) were used as inputs to the modulation classifier (block 3 of 
25 Figure 1 ). Complex baseband signals were used, whereas the carrier firequency and phase 
errors were simulated. The firequency error was random and uniformly distributed over a 
range [-4.8 kHz to 4,8 kHz], while the initial carrior phase was unifoxmiy distributed over 
[-n to 7t]. For di^ tally modulated signals, a random delay uniformly distributed over [0 to 
TsJ was used to simulate symbol timing uncertainties (where Ts is the sampling intervsl 

30 equal to 1/4S000 seconds). 
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Simulation ol Specific Signals 

Analog Modulations 

For analog modulation sctemes^ two types of fioorce signals were simulated. The first was 
a real voice signal, band-limited to [0 to 4 HTz], The second was a simulated voice signal 
ftiat used a first-order aatoregressive process of the fomi: 

ypc] 0.95 K y[k-l] + n[k] 

where n[k] is a white Gaussian noise process. Frnth^mniore, tJii? pseudo-voice signal wa£ 
band4iimt6d fuoiu 300 to 4000 Hz. 

Pot AM signals, a constant value was added to the source signal. The modulation 
index was calculated by using the maximum amplitude value over the whole soiure signal. 
The indsK was then uniformly distributed in the interval [50% to 100%] . The total length 
of the real source signal wa$ about 120 seconds, while that for the pseudo- voice was about 
40 seconds. From these two source signals, sequences of S5.3 msec in duration were 
randomly extracted. Thus^ the observed modulation index for a sequence was equal to ox 
less than the chosen modulation index. 

For jtrequency-modulated signals^ a cumulative sum was used to approximate the integral 
of the signal source. Generic FM signals were simulated using real or pseudo-voice signals 
wiili a modulation hidex unifbnnly distributed in the interv^al [1 to 4], The bandwidth 
occupied by these signals ranged iGrom 1 6 kHz to 40 kHz^ using the fipproxunation: 

where jB is the modulation index andj^ax is the maximum source frequency (4 kHz in this 
case). The AMPS FM signals were approximated nsmg a modulation index of 3, 

Digital Modulations 

Continuous-phase FSK signals were simulated hy using filtered M-ary symbols to 
provide .frequency modulation to a carrier finequency. Pager signal parameters were ba^ed 
on observations of real signals, xvith 2FSK modulation at a bit-rate of 2400 bps, a 



29 



frequsD cy deviation of 4.8 kHz, aad almost no filteaing. 4FSIC signals were also sTTiiulated, 
using the same 4,S kHz firequency deviation and a symbol rate of 1200 baud. Tbe 19.2 
kbps, 2FSK signals from the Racai Jaguar radio were simulated, using a frequency 
deviation of 6.5 kHz and a 5^ order Buttca^^OTthpre-modulalion filter with a cutoff 
frequency of 9.6 kHz, Also included in the sixnulation were GMSK signals that were 
samilar to GSM sisals ha^dng a BT product of 0-3. 

For PSK and QAISl signals, the symbols were filtered with either a raised cosine 
function or a square root raised cosine function. The selection was randoraJy performed 
with equal probabilities. The TollofT factors of 20%, 25%, 30%, SSVa, 40%, 45%, and 50% 
were imiforinly and randondy selected. For all these signals, the symbol rates were chosen 
randomly between 16 and 20 kbaud. Also simulated were -k/A-QPSK dguals tliat were 
similar to IS-54 signals, with a symbol rate of 24 kbaud and a sqxiare-Toot raised cosine 
pulse-shaping filter with a roUoff factor of 35%, 

Additive Noise 

The simulated signals were passed thiou^ an additive white Gaussian iioise 
chaimel before being classified. For 'the simulattoii, no filtering was done at the receiver, 
therefore the signal obseort'ed by the modulation classifier was corrupted by the white 
noise. The noise power was calculated from the knowledge of the average power of the 
modulated signal and the SKR over the sampling band^adth. This SNR w^ defitied as: 

SNR ^nip = S/(No . Fs) 

w^here S is the signal power. No is the white noi$^ power spectral density, and is the 
sampling frequency equal to 4S kHz, 

For amplitude-modulated signals [AM, DSB-SC, and SSB]^ the amplitude power 
was calculated by using all source signals (real and pseudo-voice^ This condition irrrplies 
that, for a given sequence, the observed SNR might be different fi^om the overall SNK, 
which is especially tme for DSB-SC md SSB signals where some segments of the signate^ 
because of silence segments, may have no pow^e{r at aE. 
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Binaiy Decision Thresliolds 

WMi respect to tilie decisLon tree analyses perfgrmed under tiie pr&sent invention, 
each decision compares a signal feature wilii a threshold. For the simulations undertaken, 
the thresholds were set from the resrults of direct obsen^tions of the feature distributions 



5 in a tramiiig set of simalated signals havang an SNR of 5 dB. These selected thresholds are 
summarized in Table 1 for die different fealnres. 



Feature 


Threshold 


Power levri above noise level, for detection of signal presence (Step 16) 


10 dB 


Power level above noise level, for bandwidth and centroid estimation 
(Step 20) 


5dB 


Maximum of the normaUzed squared FFT of the centered normalized 
envelope (U^s,\) (Step 23) 


1.44 


Variance of the direct pha^e (rad) (CWvs, FM/PSKin Step 41) 


0.25 


Variance of the direct phsi^e (rad) (AM vs. DSB-SC/BPSK in Step 68) 


4.0 


Variance of the "absolute" phase (rad) (one-dimensional vs cwq- 
dknensienal signals in Step 67} 


0.144 


Kmtosis of the instantaneous frequency (FM vs FSK in Step 44) 


2.5 


Variance of tlie normahzed amplitude in one-dimensional signals CDSB- 
SC vs BPSK in Step 70) 


0,25 


Variance of the normalized amplitude in. A/PSK signals (OTHRR vs 
QPSK/tc/4-QPSK/jWSK in Step 73) 


0-15 



Table 1; Decision thresholds used in the simulations. 



10 CI assification Results 

Aq e^tiTuate of the perfonnance of the modulation classificatioii method of the 
present invention was obtained by applying the simulated signals de$cribed earlier. For 
each modulation type, the 500 gejierated seqa^ces were classified by using the preferred 
decision tree mediods iUustr&ted in Figure 11. Eleven outputs £ram the modulation 

15 classifier were possible: NOISE, CW, AM, DSE-SC, FM, FSK, BPSK, QPSK, 7ry'4-QPSK, 



MPSK, and OIHEK Table 2 shows the classification results for each of the siimilatsd 
modulation types at a SNR of 5 dB, The resiilts of this table were obtained in the presence 
of iinifoimly distributed fi-cqucncy errors between - 4.8 kHz and 4.8 kHz (10% of the 
^ampliag frequency). The inttial carrier phase was also iraifoxmly distributed over the 
range of-iz to For digitally-modulated sigaals, a random delay uniformly distributed 
over one sample period, was used to siuiulate symbol timing uncertainties. 
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Table 2 - ClassiScatiott results (in %) for a SNR of S dB. 

In Table 2, it is impoitaat to note the difference between the classification for 
analog modulatioii signals using real voice (denoted by "V") and those itsing simulated 
voice (denoted by "SV"). For tlTese modulations, the real voice signals included pause and 
silent durations, whereas tlie simulated voice signals were generated continuously, 
■witho-tit any silent duration. For analog modulation, the pauses in. a real voice source 
produced unmodulated seq^uences. For DSB-SC and SSB modulated signals, such pauses 
produced no signal at all. In the case where the duration of the S5.3 msec sequence was 
mostly a pause, the signal was classified as noise or as OTHER, In the case where the 
transition between a pause md voice vyas not clear within tlie 85,3 msec observation, 
erroneous modulation types appeared (as the row for SSB(V) in tabic 2 indicates). For AM 
and FM signals, The absence of a source signal produced a CW classification. If the 
duration of the silence occupied most of the S5.3 msec observation time, such quiet 
sequences were classified as CW signals. Again, depend^?ig on the tra^sitinn timn hetiwen 
a pause and a voice, strange results appeared (as the row for FM(V) of Table 2 indicates). 
These results are not considered classificatidn errors as such, Radier, they reflect the fact 
that the classification is more difficult to perform for analog-modulated signals when 
usmg a very short observation time. Such perceived errors would be ehnmiited by either 
increasing dae lengtli of the observation time before the classification is undertaken, or 



33 



adding a post-processing step that acCTm^t&s the results of several observations and 
perfoffiis a decision acc^brding to the dominani modalation type. 



Modular Construction 

The modular nature of tlie clas$ificatioii process illustrated iii Figure 1 1 allows for 
different modulation scliemes to be classified. Furthermore^ ^ilthough the disclosed 
approach imder the present iuvcation assumes an AWGN clmnnel, tiiere are possibilities of 
extending the chaonel to more complex models by employing additional processing, such 
as the method of blind equalisation. 

Changes and modifications in the specifically described embodirnents can be 
carded out without departing troxn the scope of the invention which is intended to be 
limited only by the scope of the appended clairns. 
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WHAT IS CLAIMED IS 



L A method for rccognizijig the type of modulation of a complex baseband signal 
ejdxacted from a modulated signal having a earner frequency, said method 
comprising the steps of: 

(a) gen?ratiiig a pre-processed isignal fi^om the baseband signal; 

(b) examining amplitude variation? in the pre>-processed signal to Ldenti:fy the type of 
envelope thereof as being one of a constant eavelope and an irregular onvelope; 

(c) estimating the caixier frequency error in the pre-processed signal- 
ed) cgixecting the pre-processed signal for the carrier frequency errors to g^erate a 

cairier-corrected signal; and 
(e) classifying the type of modulation in the cairier-corrected signal, based on the type 
of env&lope identified in the examining step. 

2. The method of claim 1 , wherein the generating step includes a step of detecting 
presence of the baseband signal by estLtnating the power spectral density thereof ijvith the 
aid of externally measured background noigcpowcr^ md comparing said power s-pectral 
density agamst a power threshold derived firqm the background noise and a predeiined 
probability of detection. 

3 . The method of claim 2, wherein the generating step furtlier includes, after the 
detecting step, a process of gross error correction of the carrier frequency, 

4. The method of claim 3, wherein the gross eijor correction process includes the 
steps of: 

(a) estimating baseband frequency bandwidth by comparing the power spectral density 
against a bandwidth direshold derived from the background noise and a predefined 
bandwidth estimation error and a predefined probabihty of detection^ 



(b) estimating main cenbroid firequency of tilc baseband sigaaJ, by deteraoinmg a 
cottipotient centroid frequency of at least one frequency componont in tiie power 
Spectral density bearing a power above the bandwidth, threshold; and 

(c) frequency traaslating the "baseband signal by the main centroid firequency, follow^ 
by low-pass filtering with a filter bandwidth derived firom the baseband -frequeiicy 
bandwidth, 

5. The meffaod of claim 1, wherein when the type of envelope is identiJiied as a 
constant envelope, the classifying step classifies the type of modnlation according to the 
phajse and freqnency contents of the pre-processed and carrier-corrected signals, 

6. The method of claim 1^, wherein wh^ the type of envelope is identified as an 
irregular envelopes, the classifying step classifies the type of modnlation according to the 
phase and amplitude contents of the canier-coirected signal. 

7. The method of claim 1 ? wherein when the type of envelope is identified as a constant 
envelope, the step of estimating the carrier fi-equency error includes the steps of: 

(a) obtaining the npimahzed squared amplittide of the cairier-corrected signal and 
applying thereto a Fast Fotirier Transfoon (FFT) with zero-padding to generate an 

(b) squaring the absolute vaJnes of the EFT output to generate a frequency spectrum; 

(c) searching the fi-equency spectrum to find a maxknnm-power firequency sample; 

anjd 

(d) performing a fine search far the carrier fi-equency erxoi, by applying one step of a 
secant optiixdzation process to the maximum-power frequency sample and the pre- 
processed si^al- 

S. The method of claim 1, -whereia when the type of envelope is identified as an 
irregular enyclopCv ihe step of estimating the carrier frequency error includes the steps of; 
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(a) obtaining the noimsLlized squared amplitude of tbe camer-corrected si^al and 
applying thereto a Fast Fourier Transfomi wtli zero-padding to generate an FFT output; 

(b) obtaining tlie square of tlie absolute values of the FFT output to generate a frequency 
spectrum; 

(c) searching ihe j&equency spegtnun to iind a maximum-power frequency sample; 
and 

(d) performing a fine seaxcli for the carrier flrequ&icy error, by applying one step of a 
Secant optimisation process to tlie maximum-power frequency sample and the output 
of the ncamalized squared amplitude signal- 

9 - The me^thod of claim 1 , wherein when the type of envelope is identified as a constant 
envelope, the step of classifying the type of modulation mtludes ttie steps of: 

(a) ohtaining the direct phase variance for samples of the carrier-cqixected signal above 
the mean amplitude of tlie canier-corrected signal; 

(b) compamig the direct phase vmiance to a predefined phase thneshold; and 

(c) classifying the type of modulation as Contmuous Wave when the direct phase 
variance is below the phase threshold, and as Frequency Modulation when the direct 
phase variance is above the phase thTeshold. 

10. The method of clahn 9, wherein when the type of modulation is classified as 
Frequency Modulation^ the classifying step furdier includes the steps of obtaining the 
instantaneous frequency distribution of tlie pre-processed signal and obtaining the FFT of 
the 5aid instantaneous firequeaicy distribution, such that the type of modulation is classified 
as analog FM when the presence of a tone is detected in the FFT of the instantaneous 
fequency distribution. 

1 1 - The method of claim 1 0^ wherein when no tone is detected in the FFT of the 
instantaneous Jrequency distribution, the classifyiug step farther includes the steps of 
determining the iloirtosis coefficient of the instantaneous frequency distribution, such that 
tlie type of modulation is classified as Frequency Shifi: Keying (FSK) when the kurtosis 
coefficient is below a predetermined coefficient threshold. 
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12. The metbcd of claim 1 1, wheanein when the kurtosis coefficient is above a 
predeteniiined threshold, the classifying step faxQier inclades obtaining the FFT of the 
squared value of the pre-processed sigaal and cowting the number of peaks contained in 
said FPT, such that when the number of peaks is more than one^ the type of modulatiou is 
classified as FSK. having anvunber of modulation Levels e<iaal to the number of peaks; 

and as acialog FM otherwise. 

13. The method of claim 12, wherein tlie step of obtaining the instantaneous frequei^oy 
distribution includes the steps of: 

(a.) computing the instantaneous phase distribution of the pre-processed signal; 
(b) estimating the bandwidth of the instantaneous phase distribution; 
(g) low-pass filtering the instantaneous phase distribution to generate a filtered phase 
distribuliODj imd 

(d) estimating time derivative of the filtered phase distribution. 

1 4, The method of claim 1 , wherein when the type of envelope is identified as an 
irregular envelope^ the step of classifying the type of n:iodulation includes a step of 
obtaining froin the caixier-coirected signal the variance of the absolute phase 

(i) sr ^(]/(^)] + 7[2(0l) signal samples above the. mean signal amplitude, such that 

the type of modulation is classified a? beingt 

(a) one of Amplitude Modulation (AMX Double Sideband Suppressed Carrier (DSB- 
SC), and Binaiy Phase Shift Keying (BPSK), when the absolute phase variance is 
below a predefined phase threshold; and 

(b) one of Quartemaiy Phase Shift Keymg (QPSK), -:r/4-QPSK, M-axy PSEC, and 
OTHER^ when the absolute phase variance is above said phase threshold- 

1 5, The method of claim 14^ wherein when the absolute phase variance is below the 
phase thi^feshold, the classifying step further includes a step of obtaining from the earner- 
corrected signal the direct phase variance for signal samples above tlie mean signal 
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aanplitude, such thai; the type of modulation is classified as being oae of DSB-SC and 
BPSK when the direct phase variance is above the phase threshold, and as AM otherwise, 

16. The method of cLaim 1 5, wherein when the direct phase vaii ance is above the 
5 phase threshold, the classiftdng step ftirther inoludes a step of detenninhig the amplitude 

v^iriaiice of the camer-coirected signal, siuch that the type of modulation is classified as 
DSB-SC when the ampHtude variance is above a predefined amplitude threshold, and as 
BPSK otheradse. 

1 7. The method of claim 14, wherein when th.e absolute phase variance is above the 
phase tln'eshoklj the classifying step further include^ a step of obtaining firom the pre- 
processed signal the amphtade variance^ such that the type of modulation is classified as 
OTHER when the amplitnde vaiimice is above a predefined ainplitude threshold^ nnd as 
PSK otherwise. 

1 8. The method of claim 17, wherein when the amphtude vaTiai:vce is below the 
amplitude threshold, Ihe elassify^ing step further includes a step of applying a fourth power 
non-Mneaiity to the pre-processed signal followed by compuiing Fast Fonrier Transform to 
generate an EFT output; such that the type of modulation is classified as QPSK if the EFT 
output bears only one discrete component, a$ 7r/4-QPSK when the FFT output bears two 
discrete components^ and a$ M-ary PSK with M larger than 4 otherwise. 

19. A inethod for recognizing the type of modulation of a complex baseband signal 
extracted fiom a modulated signal having a carrier frequency^ said method comprising the 
steps of: 

(a) detecting presence of the baseband signal by estimating the power spectral density 
thereof with the aid of externally measured bfljcJcground noise power^ and comparing 
said power spectral density against a power threshold derived fiom the backgronnd 
3 0 noi se and a predefined probability of detection; 
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estimating baseband frequency bandwidth by comparing the power srpectral density 
against a bandwidth threshold derived from the background noise and a predejBned 
batidwidth estimation error and a predefined prdbabihty of detection; 

estimating main centroid frequency of the baseband signal by computing a 
component centroid frequency of at least one fiequency component in the powex 
spectral density bearing a power ^ove the bandwidth threshold; 

performing gross carrier error correction by frequency translating the baseband 
signal by the main centroid frequency, followed by low-pass filtering with a filter 
bandT^ridth derived from the baseband frequency bmdwidth, thereby generating a 
pre-processed signsJ from the baseband signal; : 

examining amplitude variations in the pre-proceped signal to identify the type of 
envelope thereof as being one of a constant envelope and one of an irre^lar 

envelope: ! 

< 

estimating the earner frequency error in the pre-processed signal 
by performing the steps of: * 

(i) Obt aiding the wrm^Ii^ed squared aj^iplitude of the carder-corrected signal and 

applying thereto a Fast Fourier Transform with zero-padding to generate an FFT 
output; 

(ii) squaring flie absolute values of the FFT output to generate a frequency 

(ill) searclung the frequency spectrum to find a TnaximaiEL-power fireqaency 
sample; and 

(iv) performing a fine search for the carrier frequency error, by applying one step 
of a secant optimization process to the maximurh-power frequency sample and the 
pre-processed signal; 

wlien the type of envelope is identified as a constant envelope; and 
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by performing the steps of; 

(v) obtaming the nonnalized squared amplituide of the camer-cotrected signal and 
applying ther^o a Fast Fourier Transfonn with zero-padding to generate an FFT 

5 output; 

(vi) obtaining the square of the absoluts values of the FFT output to generate a 

(vii) searching the frequ^ency specfcojm to jSnd k maximnm-power frequency sample; 
(viii) performing a fine search for the carrier feequency error, by applying one step of 

I 

10 a secant optimization process to the maximum-power frequency sample and the 

output of the normalized squared amplitude signal, 

4 

when tlie type of envelope is identified as an irregular envelope, 

1 5 (g) correcting the pre-processed signal, for the carrier frequency errors to generate a 
carrier-coiTected signal; and 

I 

(h) classi:&dng the type of modulation in the carrier-correeted signal, based on the type 
of enyelopCj, 
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by performing the steps of: 



(i) obtaining the direct phase variance for saroplcs of the carrier-corrected signal 
above liie mean ampHtude of the carrier-corrected signal; 
25 (ii) comparing the diroct phase variance to a pradcfined phase threshold; and 

(ill) classifying the type of modulation as Contmuous Wave when the direct phase 
variance is below the phase threshold^ and as Frequency Modulation when the 
direct phase variance is above the phase threshold^, 

I 

30 when the t>rpe of envelope is identified a constant envelope, and by perfqrmmg thi> step 
of classifying the type of modulation includes a step of obtaining from the carxier- 

coirected signal the variamoe of the absolute phase = ^{[/(Ol + for signal 
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samples above the mean signal amplitude^ to classify the typo of modulattoxi whom tho 
absolute phase variaace is below a predefined ptiase^ threshold as being one of Amplitude 
Modulation^ Double Sideband Suppressed Carrier, and Binary Phase Shift Keying, and. 
when the absolute phase variance is above a pre-defined phase threshold^ ais being one of 
Quartemar^^ Phase Sbift Keying (QPSK), 7t/4-QPSK, M-aiy PSK, and OTHER, when the 
type of envelope is identified as an irregular envelope. 

20- A system fot recognizing the type of modulation of a modulated signal having a 
carrier ireqnency comprising: 

(a) a receiver section for extractiiig from the modulated signal a complex baseband 
signal; 

(b) a pre-clsssification stage for generating a pre-processed signal from the baseband 

(o) means for exsmiining amplitude variations in the pre-processed signal to identify- the 
type of envelope thereof as being one of a constant envelope and an irregular 
envelope; 

(d) means for estimating the canier firequency eriror in the pre-processed signal^ 

(e) means for correcting the pre-processed signal for the carrier frequency errors to 
generate a carrier-corrected signal; and 

(f) xneans for classifying the type of modulation in the canier-corrected signal^ based on 
the type of envelope identified in the examining step. 

21 . The system of claim 20, wherein the generating means include means for detecting 
presence of the baseband signal by estimating the power spectral density thereof with the 
aid of externally measured backgroimd noise power^ and comparing said power spectral 
density against a power threshold derived from the background noise and a predefined 
probability of detectian. 

22. The system of claim 21, wherein the geaerating means ftirther includes means for 
gross error correction of the carrier firequeaicy. 
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23 , The system of claim 22, wherein wlien the type of envelope is identified as a 
constant envelope, the type of modulation is classified according to the phase and 
frequency contents of the pre-processed and carrier-carrected signals. 

24. The system of claim 22, wherein when the type of envelope is identified as an 
irregular envelope, the type of modulation is classifibd according to the phase and 
amplitude contents of tlie carrier-coirected signal. 
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ABSTRACT 

i 
I 

Disclosed is a unique systean and metliod &>t recognizing the type of modulalioii 

4 

embedded in an unknown complex baseband sigikl, comprising a receivex seption for 
extracting the complex baseband signal fiom a xnodiilated signal having a carrier 
frequency^ and comprising an orderly series of signal pK^cessing ftinctio^is for (a) 
estimating the bandwidth of the unknown signal, (b) removing the out-of-band Tioise and 

! 

correcting gross carrier frequency errors, (c) discriminatmg between constant envelope 
and irregular envelope signals, (d) estimating and correciing residual carrier frequency 
errors, (e) classifying a constant envelope signal into one of the foUomng modulation 
formats: {Continuous Wave (CW), Frequency ModiLbtion (FM), Frequency Shift Keying 
(FSK)}, a-od {£) classifying an irregular envelope signal into one of the following 
modulation format?: {Amplitude Modulation (AM), Doable Sideband Suppressed Carrier 

I 

(DSB-SC), Binary Shift Keying (BPSK), Quaternary Phase Shift Keying (QPSK), 7c/4' 
shifted QPSE. M-aiy PSK (MPSK), and OTHER classes}. 

I 

I 



I 
I 
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